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We calculate the effects of the target mass on deep-inelastic electron scattering on targets of spin 
one such as the deuteron, focusing on the novel structure functions that enter. It is important to 
understand these effects given the recent Hermes measurements of the bi{x, Q^) structure function 
of the deuteron. We also derive the spin-one analogues of the Wandzura-Wilczek relations and 
discuss possible calculations in lattice QCD. 
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I. INTRODUCTION 

' The deuteron is often used as a target in deep inelastic scattering (DIS) experiments since, neglecting nuclear effects, 
, its Fici^2d{x, Q^) and gid,2dix, Q^) structure functions are simply related to those of the proton and neutron. However, 
' the deuteron is more interesting than this; being spin-one, it has additional structure functions that cannot be reduced 
_ ', to those of the constituent nucleons W\. Recently, the Hermes collaboration used a tensor polarised deuteron gas 
target to obtain the first measurements of bid{x,Q^) in the range 0.01 < a:; < 0.45 and 0.5 GeV^ < < 5 GeV^ 
0- Their results show it to be significantly different from zero, rising as x — > 0, and consistent with coherent double 
scattering model predictions 0,01^. However, these measurements were necessarily at a relatively low and power 
corrections may be important in the interpretation of the data. 

Corrections that scale as powers of MVQ^ (where M is the target mass) are known as target mass corrections 
(TMCs) and many years ago Nachtmann |y| showed that it was possible to take these into account exactly in unpo- 
larised electron-proton DIS using the representation theory of the Lorentz group. Georgi and Politzer 0,13] extended 
I ' and clarified this, also discussing the effects of quark masses. More recently, target mass (TM) effects have been 
\ investigated in the polarized nucleon | ^ liol fill ll^ and parity violating [sl fO. IT^ structure functions, in nucleon 
■ generalised parton distributions [l^ llSl lla. Il7j and in the virtual photon structure functions j^^- Other power 
corrections arise from dynamical higher-twist effects and are much more difficult to analyse (even their separation 
from QCD scaling violations is ambiguous [3). Here we ignore these latter terms and explore the TMCs for the 
leading twist structure functions of the deuteron and other spin-one targets. We derive spin-one analogues of the 
Wandzura-Wilczek relations between moments of various spin-one structure functions and also discuss how recent 
developments in lattice QCD will enable information on the underlying parton distribution functions (PDFs) in the 
deuteron to be extracted directly from QCD. 
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II. DEEP-INELASTIC SCATTERING ON SPIN-ONE TARGETS 

The formalism of DIS on spin-one targets was set out by Hoodbhoy, Jaffe and Manohar 1] whose nomenclature we 



' attempt to follow. We define the parity conserving hadronic DIS tensor and structure functions through 

W'^''{p,q,E) ^ i- Jd'xe^^-^p,E\[J^{x),rm\p,E) (1) 

+ -5ie^"^^'^<ZAS. + ^g2e^'''^qx [p ■ qs. - s ■ qp,) , (2) 

where is the electromagnetic current and we have ignored terms involving q^ or since the DIS leptonic tensor, 
Ifiu, with which W^^ is contracted satisfies q^t^i, = q'^t^n, = 0. In Eq. |2Il the various tensors are defined as 

r^"" = ^ (^q-E*q-E-^,,^K^g'^\ 
= ^{q-E q-E--v t 



f"" = -l-fq-E* [pf^E" + p'^E''] +q-E [p^'E*" + p'' E*''] ~ -vp'^p" 
2v^ \ 3 
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the target polarisation is i?'' normalised such that — —NP and the generalised Pauli-Lubanski spin vector is 

1 + ^^,^2^ where v = p ■ q and x = /2p ■ q is the Bjorken variable (Q^ 



Finally k = 1 



andp - E ^ p - E* = 0. 

In addition to the and gi.2 structure functions familiar from spin-half targets, four additional structure func- 
tions, 61,2,3,4(2;, Q^), appear. In the parton model, 61,2(2;) have simple interpretations; defining q™^^{x) as the prob- 
ability of finding a parton with spin up (down) along the beam axis in a target with spin m along the beam axis, 
61 (x) = b2{x)/2x — ^[2q^{x) ~ q^{x) — (/{(s;)] 0- Since the accompanying tensors are symmetric in their indices 
and vanish when averaged over target polarisation ((-E^-Ej^)spin avg. = \—M'^g^^w /3), these functions are mea- 

surable with unpolarised lepton beams but require a tensor polarised target. Various models of the novel structure 



functions have been proposed [J, U, |5|, |21|, |22, 123, |2J| and bounds arising from positivity are known j25|, |2y| . 

Using the optical theorem, the hadronic tensor above is related to the forward Compton scattering amplitude T^'', 
= ^%Ti"']. We decompose the tensor structure of T^'^ as for 1^'^'' in Eq. O with {Fi,g,,b,}. At 

leading twist, the spin-one Compton amplitude is given by {ipf{x) is a quark field of flavour / and electric charge Qf) 

T^-"- ^{p,E\t'^''\p,E) = I j d^xe''^--'{p,E\TY,Q}i^fri'f{^)i^fl''i'fmP,E) + ^ J^, q ^ ~q] (3) 



71^2 

even 



9^1 



n-l 



+2zE>V(3)e 



i9A 



n=l 
odd 



2 



q^Jp,E\0^^-^-\p,E) 

q^Jp,E\0^'^-^^"\p,E) 
q^ Jp,E\d''^ \p,E), 



(4) 



where the wi*'' (Q^, /i^, Ofg) are Wilson coefficients (which are identical to those in spin-half DIS), fi is the renormali- 
sation/factorisation scale, and the vector and axial- vector twist-two operators are 



(5) 
(6) 



Here D^^ = ^{D — DY, {. . .} indicates symmetrisation over the enclosed indices and traces are subtracted such that 
the operators transform irreducibly under the Lorentz group - they have fixed spin. The matrix elements of these 
operators in a spin-one target of momentum p and polarisation E are given by 



{p,E\0^'^■■■^'-\p,E) = a„ [p''^ . . .p^" - tr] + d„ 



3^ ^ 



tr 



(p,i?|0''-''"|p,i?) = 



tr 



(7) 
(8) 



The parameters a„, d„ and r„ (we have chosen a slightly different normalisation for r„ than Hoodbhoy et al. 
such that it is dimensionless) depend only on the choice of target and the renormalisation scale. They encode all the 
non-perturbative information about the structure of the target. The novel effects in spin-one targets arise from the 
additional structure proportional to rf„ in Eq. Q; this vanishes when averaged over polarisations even in the presence 
of the trace subtractions. 
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III. TARGET MASS EFFECTS 



Away from the Bjorken limit, the traces in Eqs. and ((Sj) can not be ignored, but by including them we can 
treat TMCs exactly 0, ■ The contractions of the trace subtracted tensors required in the Compton amplitude are 
given explicitly in the Appendix, Eqs. ljA3|l - ljA7|l . and inserting these results in Eq. Q leads to target mass corrected 
expressions for the various functions {Fi,gi,bi}. We find 



where 



F,{p,q) ^ Y: 2a„(^) |w«C«(^)--^-^W(^)c(^2,(r;)^, (9) 



n— 2,4,... ^ / \ / 



Up,,)- - E ^^^(I) + + ^ (12) 



n=l,3,. 



(10) 



Mp,,)= E ^^^(|)"{^^i^-^i(^) + (- + 3),ci^^^ (11) 



^3(p,«) - E (^)"" {8'y^cil('y) + Vcit^3('?) - , (15) 



^ VP <1 

and the various Cn "* (??) are Gegenbauer polynomials |23] ■ 

The TMCs in these formulae can be removed by taking only the leading power of ry in the curly braces to give 
inverse powers of x = Q'^ /Aprj; in this limit the above forms reduce to those given in Ref. Q. The absence of terms 
proportional to dn in Fi,2{p, q) (which results from non-trivial cancellations) could be expected since these correspond 
to the polarisation averaged structure functions. Also, the Callan-Gross relation, F2{x,Q^) = 2xFi{x,Q^), and its 
spin-one analogue, b2{x, Q^) = 2xbi{x, Q^) are broken by target mass effects even at leading order in QCD where 
Wk'^ = 1. 

To determine the target mass corrected structure functions, we follow the method of Georgi and Politzer 0. We 
will outline this explicitly for bi and give the final results for the other structure functions. Using the series expansion 
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of the various Gegenbauer polynomials [23 changing the orders of the subsequent summations, we can rewrite 
Eq. lO as 



£=0 j=0 



''^21+2] +2^21+2 j +2 



AP M%2j+4d2i+2j+4 [21 + J + 3) 

(2^ + 2j + 2)(2^ + 2j + l) 02 (2^ + 2j + 4)2(2£ + 2j + 3)2 



(2€ + j + 2)! /M^y 



(18) 



where uj = 1/x and we have suppressed the and dependence of &i for brevity. The coefficient of oj" ^ in this 
expansion enables us to extract the Mellin moments of via 



MJbi) = [ dxx"-'^bi(x) = — I dujuj-"bi(uj) 
Jo 2Tn J 



(19) 



(n-\ 



(1) 

n-l+2j' 



3=0 



(n-l + 2j)(n-2 + 2j) j!(n-3)! V Q' 



To proceed, we introduce the functions i^j*^ (y) whose moments are given by 



d„-i+2, (n-l+j)!/M2y w(^_^i+2,d„+i+2, (n + j)! /M^y 



i!^ +4— y^Ili 
Q2 / ^ Q2 ^ (n + 1 



Q\t^ + 1 + + 2jT - 3)! V 



whereby it is trivial to show that 



and 



^^^^= j\yy-^Kf{y)^ Where Kf{y)^fdzj\z'Hli\z'). 



Y)= dvv''-^Hf {y) , where (y) = dz j' dz' F^p {z') , 



(20) 



(21) 



(22) 



The f'P'^ and ^ ete. differ in a known, perturbative QCD manner determined by the Wilson coefficients. For use 
below, we also define 



W«a„= / dyy-F^\y), F«(y)= / dz / dz'F«(z'), 



(23) 



and 



W«r„ = / dyy-F^^ (y) , H^^ (y) f' ^ f' dz' F^^ [z') . 



(24) 



from the coefficients a„ and r„. Now noting that 



E 

J=0 



(n + j)!^j- ^ 1 
j\n\ ^ (1 - z)"+i 



an inverse Mellin transform of Eq. p9|l gives the target mass corrected structure function: 



(25) 



hi{x,Q\M') 



2Tii 



dn a;-"+^M„(6i) 



— 'too 
■2 /-l ^2 



(26) 



g2 7o "y(l - M2y2/Q2) j 27ri 
Using the identity dn z" = 27r« (5(ln(z)), this finally leads to 



dn 



y/x 



1-M2y2/Q2 



bi{x,Q^,M^) = 2a:— 



hPx^d^ 

02 



(27) 
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where 



2x 



1 + ^1 + 4Af2a;2/Q2 



(28) 



is the Nachtmann variable |a] and r — yyi + AM^x^jQ^ . 
For the other structure functions wc similarly find 



(29) 



(30) 



(31) 



g2(x,g^^f2 



'2a;—— X 
ax'^ 



(32) 



b2{x, 



4^1 2 — 



(33) 



b3ix,Q^,M^ 



dx^' 



4— X 



dx^ 



ifr^(e) 



(34) 



64 (x, ( 



dx4' 



dx^ 



dx^ 



(0 



(35) 



Since terms proportional to dn vanish identically in the Fi structure functions, the form of the TMCs in these (and the 
gi) structure functions is the same as those in the spin-half case. However, the parameters a„, d„ and r„ describing the 

matrix elements in Eqs. Q and (jSJ depend on the specific target. Setting Wn'^ = 1, our results for i^i,2(a;, Q^, M^) are 
equivalent to those of Georgi and Politzer 8] and those for 51,2(2;, Q^, AP) agree with Refs. jlltil^] when twist-three 
operators are ignored. Since the only modifications occur because of the additional structure in the matrix element 
in Eq. Q that is not present for spin-zero or spin-half targets, we can immediately conclude that the target mass 
effects in Fi and gi are of the same form in targets of any spin greater than one-half. It is also clear that TMCs in 
spin-zero targets such as the pion have the same form as those in the unpolarised structure functions of the proton. 

Since we have only kept twist-two operators in the expansion of Eq. (PJ, Mellin moments of gi and 32 above are 
connected by the Wandzura-Wilczek (WW) relations even after target mass corrections are included 0, ^3 • 
Similarly, our results for 62.3,4 satisfy the relations: 



1 71 — n — 6 



— n 



M„(&2), 



M„{bi) = 2 M„(62), 

3 — n 



(36) 
(37) 



where M„(/) = c?xx"~2/(x). These relations are also valid in the presence of target mass and perturbative QCD 
effects and can generally be expected to hold to better accuracy than the usual WW relation since, for massless 
quarks, additional operators only contribute at twist-four (as opposed to twist-three in the WW case). 
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FIG. 1: Structure functions of the proton. In the main panels we show the F2p (left) and gip (right) structure functions of the 
proton for Mp = 0.938 GeV (solid curve) and for Mp — (dashed curve) at a scale = 3 GeV^. In the insets, the ratio of 
structure functions for Mp = 0.938 GeV to those for Mp = is shown for — 1, 3, 10 GeV^ (dotted, solid and dashed curves, 
respectively). 



IV. DISCUSSION 



The results presented m the preceding section are relevant for the data on bid{x,Q'^) obtained by the Hermes 
collaboration since they work at a relatively low momentum transfer where 0.75 < Mj/Q^ < 6.9. If the data is 
analysed to extract the twist-two parton distributions in the deuteron {q^{x), ql{x) and q^{x) in the notation of 
Hoodbhoy et al. 01) using the parton model result, bi{x) — ^[2q^{x) — — g|(a;)], or a QCD improved form, these 
TM effects will lead to significant modifications. In order to assess the impact of the TMCs, we first compare the 
effects in F2{x, Q^) and gi{x, Q^) in the proton and the deuteron and then discuss the new results for 6id(a;, Q^)- 

Figure n shows the F2p and gip structure functions of the proton at the physical mass, Mp = 0.938 GeV, and in 
the limit Mp —> 0. To explore the TM effects in a simple manner, we work at zeroth order in the QCD coupling. 

As input, we have set Wn^ — 1 and used the unpolarised parton distributions of Alckhin 28] and the polarised 

PDFs of Bliimlein and Bottcher (Scenario 1), both at leading order: Fp\x) — [4:u{x) + d{x) + 5q{x)] and 

f!:^\x) = [4AM(a;) + Ad{x) + 5Aq{x)]. In the left panel of the figure we show F2p{x, ^ 3 GeV^) from Eq. ^ 
evaluated with Mp = 0.938 GeV (solid curve) and Mp — (dashed curve). In the inset we show the ratio of these 
functions for =1, 3, 10 GeV^ (dotted, sohd and dashed curves, respectively) . In the right panel, the same curves 
are shown for the polarised structure function gip (see Ref. [ll|, ISfl for similar studies) . One sees the well known 
phenomenon that as x — > 1, the target mass corrected functions do not vanish as they should kinematically. However 
higher twist effects that we have ignored in the operator product expansion of Eq. become increasingly important 
as X approaches unity and considering only TMCs in the large x region is not reliable 0, |^ l3^ l33l l33 | . 

In Figure 12 we show similar results for the deuteron structure functions i^2d and gid, comparing the structure 
functions for Md = 1.876 GeV to those in the zero target-mass limit. Again we use the PDF parameterisations of 
Refs. mmi, setting fP{x) = [u{x) + d(a;) + 2q{x)] and F^^\x) = ^ [Au{x) + Ad{x) + 2Aq{x)]. The larger 
mass of the target and the different flavour composition of the deuteron leads to slightly larger target mass effects in 
this case. In all figures, the abscissae are defined with respect to the proton mass, Xp = /2MpqQ. 

It is less clear how to quantify the TM effects in the case of the structure function since we do not know the 
underlying twist-two parton distributions from data at large . As an estimate, we choose two forms that are fit to 



Here we have ignored small nuclear effects in the deuteron, assuming it to be a free proton and neutron. However for heavier targets 
of spin one and higher it is interesting to consider the analogue of the EMC effect. For the 61 structure function, one may look for 
the deviation of the ratio feiyi(x)/bi(j(a:) from unity for a series of nuclei of spin-one, e.g. ®Li, ^''N. This will be non-zero and will be 
unrelated to the EMC effect observed in F2 structure functions. 
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FIG. 2: F2d (left) and gid (right) structure functions of the deuteron. In the main figures we show the Md ~ 1.876 GeV (solid 
curve) and Md — > (dashed curve) results at = 3 GeV^ . In the insets the ratio of physical- to zero- mass structure functions 
is shown for = 1, 3, 10 GeV^ (dotted, solid and dashed curves, respectively). 



Xpbid 0.001 




FIG. 3; Calculated bid structure functions of the deuteron. In the left plot we use the parameterisation of Eq. 13811 whilst 
the right plot uses Eq. (gnj- In the main figure we show Xpbid for Md = 1.876 GeV (solid) and for Md (dashed) at 
= 3 GeV^ along with data from the Hermes collaboration (diamonds). In the insets the ratio of physical- to zero- mass 
functions is shown for — 1, 3, 10 GeV^ (dotted, solid and dashed curves, respectively). We also display the relative shift 
at the experimental kinematics (diamonds). The kinks in the inset of the right plot are an artefact of parameterisation (II) 
having a zero that is shifted by the TM corrections. 



the data using minimisation: 



(I) 
(11) 



Aix'"{l-x)'^ 



(38) 
(39) 



where the fit parameters are Aj = 0.00159, 6/ = -1.828 for fit (I) and An = 0.0586, bn = -1.532, c// = 2.86, 
dii ~ —4.32 for fit (II). These give a reasonable description of the data (see dashed curves in Figure EJ- Figure 01 
shows the structure functions that result from Eq. (|38|l (left) and Eq. H39|) (right) for the physical deuteron mass 
(solid curve) and for Md (dashed curve), both at Q'^ = 3 GeV^. In the main figures we also show the Hermes 
data (diamonds, error bars include statistical and systematic errors), and in the inset we show the physical- to zero- 
mass bid ratio at — 1, 3, 10 GeV^ (dotted, solid and dashed curves) along with the ratio evaluated for the different 
Hermes kinematics (diamonds) 0. 

From Figure 13 we see that relative shift from TMCs can be significant (especially if the twist-two PDFs is negative 
in some range of x), though in the region where the structure function itself is small. If we were to go further and 
attempt to extract the three independent iso-scalar twist-two PDFs in the deuteron from data on the three structure 
functions F2d, gid and bid, it would be necessary to include these effects. However, because of the low scales at which 
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the measurements of bid have been made, QCD effects (which are independent of the target) and higher-twist power 
corrections^ in the data will also be important. Without very extensive and accurate studies that would quantify such 
contributions, a reliable extraction is difficult. An approach that may help in this regard is to use lattice QCD to 
calculate moments of the twist-two piece of bid{x, Q^) since matrix elements of the twist-two operators, Eqs. Q and 
can be studied directly on the lattice. Indeed the various structure functions of the p- meson including bip have 
been explored [s^- For ^ two hadron system such as the deuteron, calculations are more involved. Recently it has 
been shown that calculations of the lattice volume dependence of two-nucleon energy levels in background electroweak 
fields will lead to determinations of electroweak properties of the deuteron |3^ . Similar methods, employing external 
twist- two fields, will also enable the extraction of the matrix elements in Eqs.© and © [ssj . 

To summarise, we have discussed target mass corrections in the structure functions of spin-one targets. These 
effects are relevant for the recent Hermes measurements of Q^) in the deuteron and should be included in any 

extraction of the underlying parton distributions. We have also derived non-trivial relations between the moments of 
the spin-one structure functions 62,3,4(0:, Q^) that are the analogues of the well-known Wandzura-Wilczek relations. 
The calculations described here can easily be extended to targets of arbitrary spi n and to parity violating structure 
functions. The formalism for DIS in these cases is detailed in Refs. ^39, 40] and jl^] respectively, and the calculational 
procedures are the same as those used here. However with no imminent prospect of experimental tests we leave this 
to the future. 
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APPENDIX A: TRACELESS TENSORS AND CONTRACTIONS 

be a completely symmetric tensor constructed of n vectors a^' . Then the traceless piece 

can be constructed as 



(Al) 



where the harmonic polynomial |42 



fix) = y: 



k=0 



{n — k)\ f x' 
i\k\ 



This construction enables us to calculate the contractions of the matrix elements required in Eq. Q : 

-tr] = p"CW(77), 



^*{Mi^A'2pM3 ^ ^ pM,.} _tr 



n{n — 1) 



(A2) 

(A3) 
(A4) 



s^^ip'^2...p''"> -tr 



n-l 



(A5) 



{(n + 2)i^&2^{^) + Ar^^c'i\{i^)] 



9mi • • ■ (li^^ 



■•■P""^ - tr 



2^" 



(n-f l)(n-|-2) 



(A6) 



^ A model calculation of the twist-four contributions to the lowest moment of b\ finds a small effect l35l . 



9/^1 •• • %r. 



2p" 



[(n + l)(n + 2)]2 
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(A7) 



P 
~q' 



{2Aa{q,E)Cl^lM + \Ef 



+2^{2Aa{q,E)d:iM + \E\ 
where a{q, E) — [q ■ Eq ■ E*] jq^ and in the last two relations we have made repeated use of the derivative identity 

■ ■ ■ q^n+^Ti^'-'-^'^ =q,.--- q.^ri'^"-^ . (as) 
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